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Abstract 



A quantum random walk on the integers exhibits pseudo memory 
effects, in that its probabihty distribution after N steps is determined 
by reshuffling the first N distributions that arise in a classical random 
^ ■ walk with the same initial distribution. In a classical walk, entropy in- 

^ , crease can be regarded as a consequence of the majorization ordering 

of successive distributions. The Lorenz curves of successive distribu- 
tions for a symmetric quantum walk reveal no majorization ordering 
in general. Nevertheless, entropy can increase, and computer experi- 
r> I ments show that it does so on average. Varying the stages at which 

■ the quantum coin system is traced out leads to new quantum walks, 

including a symmetric walk for which majorization ordering is valid 
but the spreading rate exceeds that of the usual symmetric quantum 
walk. 



During the stochastic evolution of a classical random walk (CRW), corre- 
lations are established among the states of its two constituent parts, a coin 
and a walker. Because of the widespread use of the CRW in applications 
involving classical computer simulations [Tj, recent interest in quantum com- 
putation has focussed attention on the notion of a quantum random walk 
(QRW). The key idea is to replace the classical correlations between coin and 
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walker states in a CRW by the emblematic notion of quantum correlation, 
i.e. entanglement |2] of the states of suitable quantum analogues of the coin 
and walker. Of main interest in QRW studies has been the effect of entan- 
glement on various asymptotics, on spreading properties, and on hitting and 
mixing times. From the earliest formulations of QRWs [Sj |1] , to recent stud- 
ies on general graphs 0, on the line [6j etc., it has emerged that a number 
of surprising features distinguish quantum from classical walks, such as their 
non-Gaussian asymptotics, a quadratic speed up in spreading rate on the line 

an exponentially faster hitting time in hypercubes jH] [HI, and an exponen- 
tially faster penetration time of decision trees JDIPI]- These findings have 
recently prompted proposals for physical implementation of such processes 
in experiments, e.g. in ion traps ^21, optical lattices or in cavity QED 
|14j . A review providing a comprehensive introduction and other references 
has recently appeared [TH] . 

There are two different ways to consider a QRW as a process. In the 
first, a fixed number N of applications of a unitary quantum evolution to 
a combined coin-walker system is considered, producing a highly entangled 
state, and after the last step, the coin degrees of freedom are traced out to 
determine a reduced density matrix and subsequently Pq^\ the associated 
probability distribution (pd) of location probabilities for the walker. In the 
second, extended interpretation, the QRW is considered as the process that 
produces successively in this way, the pds -Pq^^ for = 0, 1, 2, 3, . . . . The 
purpose of the present note is to indicate some further remarkable properties 
that differentiate a QRW, interpreted in this extended way, from a CRW. 

We construct the QRW on H = He ® H^j, where the coin space is He 
= hiiO, 1}) and the walker space is H^j = h{Z), and we consider the unitary 
evolution operator V = P+U ® + P^U ® E_ acting on H. Here U is 
2x2 unitary, and acts on H^ along with P+ and P_, which are orthogonal 
projectors onto the coin states |0) = (1 0)^ and |1) = (0 1)^ respectively. 
The step operators act on \m) G if^ for m = 0, ±1, ±2, ... as 

E±\m) = |m ± 1), so that E±Ezfi = 1. Together with the distance operator 
L acting as L\m) = m\m), they satisfy the commutation relations [L, = 
±£'-1-, [E^,E_.] = 0. With pi"-* = \^){ip\ chosen as the initial density matrix 
of the coin and pw^ = |0)(0| the initial density matrix of the walker, the total 
initial density matrix p*^*^) is 

p^°^ = pi°^®pL°^ = l^)(^l®|o)(o|. (1) 

After A^ successive applications of the unitary evolution defined by V, the 
coin degrees of freedom are traced out, so that evolves to pl!^\ given by 
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the action of a trace preserving completely positive (CP) map EyN JHI 



(2) 



A:=0,1 



k 



0, 1 by Al 
- 1. 



(TV) 



{k\V^\ip), and 
(3) 



Here the Kraus operators are given for k 
satisfy 

fc=0,l 

The diagonal element {k\p[!^^\k) of the walker density matrix, for k E {—N, —N+ 
2, . . . A^}, is the probability PQ^\k) of occupation of the site k by the walker 
after the Nth step. It follows from (P) that Pg ^(A;) = 5^0- 

/ a^N) g{N) 

With V defined as above we have V = ( ^^jy) ^{n) 

operators a^'^^ etc. acting on H^. For the general choice \ip) = c|0) + d\l), 
with c, d complex, we have that A^^ = ca^^'^ + d(3^^\ Af^ = c-f^^'^ + d6^^\ 



for suitable 



We choose U = U{p) 



, with < p < 1, without any 



significant loss of generality ||5j. The choice c = 1/v^, d = and p = 1/2 
is known [HI E] to result in P^^\-k) = P^^\k) for all k, but a QRW is 
not symmetric in general. In what follows, we utilize the elementwise or 
Hadamard product A = A o B, defined between matrices A,B of the same 
size by {AoB)ij = AijBij, and we call doubly stochastic [Ej a square matrix 
A (of finite size, or with columns and rows of finite support, as is the case 
here), with nonnegative elements, which has unit column and row sums. 
Our first result is summarized as: 



Proposition 1. There exists a doubly stochastic matrix Ag'^ = A)^'' ' o 



x(N) 



Aq'' + A^^^ o A^-^\ that connects the initial pd Pq' with the pd Pq 



(0) 



f of 

the Nth step. These matrices satisfy the inhomogeneous recurrence relation 

= ((i-p)E+ + p£;_))Af) 



+ (2^ -1)4"^ 04-)] 

AN) 



(4) 



By means of this relation, Pq" ' is related to the classical pds {P^\ Pj^ ^\.. 
Pq^}, arising from the first N steps of a CRW, via the map 5*^^^ given 



P 



(1) p(0) 



C ' 
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by 



+... + ..(^-i)p«+.;(^)p(°\ (5) 

where we have introduced the reshujfling matrices ti;(*+^) = [E^—E^)M^'^\ i > 

with M(^) = ^p{l-p){A^^ o + o 4'^) + (2p - 1)4'^ ° 4'^ 
TOt/i M(o) = 0, cj(i) = 0. 

The proof is straightforward. 

Comments: 

i omce the are polynomials of degree in the commuting step 

operators they are normal operators i.e. A\!^^ ^i^^^ = ^i^^^ ^i^^ f*^^ 
k = 0, 1. This means that in addition to Q, there is a similar relation with 
A^j!^\ interchanged. Together, these two relations lead to the dou- 

ble stochasticity of Aq^'' jTHj. This simple method of constructing doubly 
stochastic matrices by convex sums of Hadamard products of Kraus genera- 
tors of CP maps, being normal operators, together with the question of the 
ensuing entropy increase, to be discussed shortly, is a nontrivial extension of 
Uhlmann's theory which addresses those questions for unitary CP maps only 
(of. and references therein). 

2) The general form of © is A^f"^^^ = AcA^f^ + {E+ - E^)M^^\ and 
we see that the final A^-dependent inhomogeneous term, which must have 
zero column and row sums, distinguishes a QRW from a CRW and moreover 
carries the burden of possible breaking of the majorization ordering, as will 
be seen shortly. 

3) Proposition 1 shows that the effect of tracing out the coin system 
after applications of V to p^^\ results in a kind of pseudo memory ef- 
fect (or pseudo nonMarkovian effect), in that determination of the quantum 
occupation probabilities at step involves the occupation probabilities of 
an A^-step CRW, reshuffled from step to step as in (0). This suggests a 
modified QRW where the coin system is traced out after every m steps, for 
some fixed m, rather than after 1 or 2 or . . . A^ steps as in the QRW as 
considered to date. The case m = 1 defines a scheme that promptly traces 
the coin system after each V action. With U{p) as above and p*^"-* as in 
(HI), this yields the occupation probabilities of a (p, 1 — p) biased CRW. In- 
deed if pl^+^) = (ev)^+i(pS')) = SvipL''^) ^ Tr,[V{p^^'> ® pL''^)V% then 
we have at each step, along the diagonal of the reduced density matrix, the 
probabilities of the corresponding row of the classical Pascal triangle i.e. 
{k\p^w \k) = {PQ^^){k) = {Pl^^){k). In this case no memory effects are 
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present and -Pq^^^'' = AcPq^\ where Ac = (1 — p)E+ + pE_ is a doubly 
stochastic matrix that repeatedly mixes the evolving probability distribution 
and extends its support by one unit to the left and right at each time step 
of the walk. 

4) Proposition 2 below is devoted to the case m = 2, where a pseudo 
memory effect is also exhibited, since here also the determination of the 
quantum pd at step N involves classical occupation probabilities from the 
first steps. However now the reshuffling matrix $ is fixed (as in © below). 

Proposition 2. There exists a doubly stochastic matrix Aq = BqoBq+Bio 

El that connects the Nth step pd Pq^^ identified with the diagonal elements 
of p[!^\ with Pq^'^^^ at the {N + 1) th step, identified with the diagonal 
elements of pi!i^^^^ = ev^{pw'')-Here the Kraus generators are 



Bo = {pc + ^p{l-p)d)El + {{l - p)c - ^p{l-p)d)l, 
El = {pd- ^p{l-p)c)El + {{l-p)d+ ^p{l -p)c)l. 

The recurrence relation satisfied by Pq^\ and its solution, are given by 

P^' = A«i^"-' = (Ac + and 

p<;" = Ajp™' = (Ac + i>rpl," 

A:=0 ^ ^ k=0 



where Aq = Ac + $, and where the matrix $ with null column and row sums 
is given explicitly by 

$ = 5o o 5o + fii o fii - Ac 



= \pc + ^/pil-p)d\^El + I -p)c - pd\^E^_ 

+ (1(1 -p)c - ^p{l-p)d\' + 1(1 - p)d 

+ Vp(1-p)c|')1 - (1 - P)E+ - pE_. (6) 

Once it is seen that in this case the analogue of © reads p^~^^^ = Sy'^ipw '^) = 
Trc[1/2Ar^^(o)^^(o))(y2t)7V] _ ^^^^^^ E^p^^^ E^ the rest of the proof is similar 
to that for the Proposition 1. Generalization to the case m > 2 is straight- 
forward. 

It is known that CRW pds become more entropic as increases. This 
can be attributed to the fact that they are ordered by majorization pUj . 
Thus, for two consecutive classical pds P'^'' and P'^^^\ each with elements 
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arranged in nondecreasing order, it is true that >- P^^^^"*, and therefore 
that S{Pc^^) < S{Pc'^~^^^), for the respective Shannon entropies, defined as 
S{P) = — J2k -P(^) log-P(A;). To facihtate a comparison with the correspond- 
ing behaviour in the symmetric QRW of Proposition 1, we set the classical 
pd P^"* in the upper horizontal line in Fig. ^ and the quantum pd Pq^^ 
in the lower horizontal line. We find the remarkable result that in certain 
cases, though QRW pds are becoming more entropic in the course of time, 
namely S{Pq^^) < S{Pq^^^^), majorization breaks down i.e. Pq^^ )f- Pq^^^\ 
even in the early stages. This is illustrated in Fig. |21 which refers to steps 

= 6, 7, 8 and 9 of the symmetric QRW. Here the entropy is increasing, 
with S^{Q) ^ 1.6551, S^{7) ^ 1.8138, S^{8) ^ 1.8909, S^{9) ^ 1.9295, but 
majorizarion ordering is violated, as can be seen at once from the correspond- 
ing Lorenz curves. The Lorenz curve of a pd P whose elements P{k) have 
been arranged in nonincreasing order is the plot of the points {n/N, '-)n{P)) 
forn = 0, 1, . . . iV, where 7„(P) = ELi and 7o(P) = 0, 7^(P) = 1. 

If P -< P' then the Lorenz curve of P' always lies below that of P and never 
crosses it. Thus crossing of Lorenz curves implies majorization breakdown 
and vice versa. Such a breakdown of majorization takes place in the sym- 
metric QRW as is seen in Fig. 2, while entropy increases as seen in Fig. 3. 
Remarkably, this is not the case for the QRW of Proposition 2, since in this 
case the mixing matrix from step to step remains iV-independent, fixed and 
doubly stochastic, guaranteeing majorization ordering [20} l2T]. 

Fig. El illustrates three notable features of entropy dynamics common 
to three symmetric QRWs as in Proposition 1, with U{p) having p = 1/3 = 
cos^(7r/6), p = 1/2 = cos^(7r/4) and p = 3/4 = cos^(7r/3), respectively: 
firstly, an increase of entropy on average, e.g. in the case p = 1/4, the 
sequence of steps 45, 51, 57, 63, ... has monotonically increasing entropy 
values; secondly, clusters of steps with decreasing entropy, e.g. S{Pq^^) ~ 

3.3498, S{P^^^) ^ 3.3467, S{P^^^ ^ 3.3408; and thirdly, a larger rate of 
increase of quantum entropy cf. classical entropy on average. 

Finally, we indicate that a symmetric QRW with a constantly delayed 
tracing scheme, as in Proposition 2, shows an even greater rate of spreading 
than the usual symmetric QRW, which in turn is known to spread quadrat- 
ically faster than a CRW 013]. 

Define the mth order statistical moment of the the distance operator L at 
step by {L^)n = Tr{pi!^^ L^), where for the "classical" case correspond- 
ing to the "promptly traced" QRW, pU^^ = {ev)^{pw^), and for the cases 
of quantum walks described in Propositions 1 and 2, pif-* = eyiv(p^^) and 
pw^^ = {£v^)^{Pw^)^ respectively. Consider symmetric walks in each case, 
with p = 1/2. All first moments are zero in each case, i.e. {L)n = 0, for 
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all A^, so that the standard deviation at step N is given by = {l?) at. 
For the CRW we have cx^ = ^J~N as is well known pp. Our methods al- 
low us to calculate uat easily for QRWs, for a given N . For the first five 
steps of the QRW of Proposition 1 we get aj^^ = erf, cr^^ = erf, = 
0-3, crf^ = {V5/2)a^, a^^ = ^/8/5a^. The enhanced rate of growth 
in the quantum case, which is known [71 U] to be given as — > cxd by 

a}^^^ ~ xJn{2 - V2)/2a%, is soon clear. 

For the first five steps of the QRW of Proposition 2 we have 
af , af = , af = VSa^ , af = ^a^, af = 2al 

We see that the standard deviations for this 'delayed tracing' QRW grow 
even faster than those of the first type of QRW. Since this second type of 
walk has constant Kraus generators, it may well be more easily implemented 
experimentally than the first type fT^ 1131 ITl] . 
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Figure 1: Successive distributions of a CRW (upper horizontal line) obtained 
by action of Ac, and of a QRW as in Proposition 1 (lower horizontal line), 
obtained by the action of A[^\ The pseudo memory effect is shown by the 
vertical arrows 
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Lorenz Curves for Quantum Probability Distributions of Steps 6,7,8 and 9 
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Figure 2: Lorenz curves for the distributions of steps 6,7,8 and 9 of the 
symmetric QRW of Proposition 1. 




Figure 3: Quantum and classical entropies v. number of steps. The symmet- 
ric QRW is as in Proposition 1 in each case. 
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